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ABSTRACT 

In this paper, the influence of binary stars on the measured kinematics of dwarf galaxies 
is investigated. Using realistic distributions of the orbital parameters (semi-major axis, 
eccentricity, ...), analytical expressions are derived for the changes induced by the 
presence of binary stars in the measured velocity moments of low-mass galaxies (such 
as the projected velocity dispersion and the 4 th order Gauss-Hermite coefficient hi). 
It is shown that there is a noticeable change in the observed velocity dispersion if the 
intrinsic velocity dispersion of a galaxy is of the same order as the binary velocity 
dispersion. The kurtosis of the line-of-sight velocity distribution (LOSVD) is affected 
even at higher values of the intrinsic velocity dispersion. 

Moreover, the LOSVD of the binary stars (i.e. the probability of finding a star 
in a binary system with a particular projected velocity) is given in closed form, ap- 
proximating the constituent stars of all binaries to revolve on circular orbits around 
each other. With this binary LOSVD, we calculate the observed LOSVD, the latter 
quantifying the movement of stars along the line of sight caused both by the stars' 
orbits through the galaxy and by the motion of stars in binary systems. As suggested 
by the changes induced in the moments, the observed LOSVD becomes more peaked 
around zero velocity and develops broader high-velocity wings. These results are im- 
portant in interpreting kinematics derived from integrated-light spectra of low-mass 
galaxies and many of the intermediate results are useful for comparison with Monte 
Carlo simulations. 



1 INTRODUCTION 



Previous authors have used Monte Carlo simulations 
to estimate the influence of binaries on the observed veloc- 



The usual approach to derive the kinematics (mean veloc- 
ity along the line of sight, the velocity dispersion, . . . ) of 
galaxies from observed spectra is comparing them to the 
spectra of so-called template stars. The broadening of the 
absorption lines in the galaxy spectra is then interpreted as 
a result of the orbital motion of the myriad of stars along 
the line of sight through the galaxy : each star has a Doppler 
shifted spectrum and the absorption lines consequently ap- 
pear at slightly different wavelengths. These template stars 
are carefully selected : they should not rotate too rapidly, 
not be member of a binary system or be peculiar in any other 
way so as to be a good representation of the average stel- 
lar population. Thus, the broadening of the absorption lines 
is explained completely in terms of stellar orbital motions. 
However, a large fraction of the stars in a galaxy are mem- 
bers of a binary (or even multiple) system. Stars in binary 
systems orbit the center of mass of the two stars and this 
extra velocity adds to the velocity dispersion of the galaxy. 
Hence, the measured kinematics of galaxies are not indepen- 
dent of the population of binary stars. 



ity dispersion (see Hargreaves et al. ( 1996 ) and references 
therein). This work focused mainly on observations of dwarf 
galaxies in the Local Group or globular clusters that can 
be resolved into individual stars. For a sample of stars in a 
galaxy, radial velocities are measured and from these veloci- 
ties, the mean velocity and velocity dispersion profiles can be 
established. The Monte Carlo calculations by construction 
are very suited to interpret these "discrete" data, e.g. it is 
possible to check the effect of repeated observations to weed 
out short-period binaries. The crowded inner regions of glob- 
ulars or dwarf galaxies outside of the Local Group however 
are impervious to this kind of study. There, integrated-light 
spectra are needed to measure the kinematics. As bigger 
telescopes make it possible to penetrate to ever lower lev- 
els of surface brightness and hence to study low-mass stellar 
systems at large distances, it is important to investigate the 
effect of binary stars on the derived LOSVD. Since Monte 
Carlo calculations at present are unable to yield the full 
LOSVD or even higher order moments such as the kurto- 
sis without excessive computational effort, they cannot be 
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Figure 1. The orbit of the primary around the binary's center of 
mass. Indicated are the descending node DN and the ascending 
node AN, the longitude of ascending node Q with respect to some 
reference direction p, the argument of the periastron u>, the true 
anomaly <fc and the inclination i. 
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Figure 2. The specific binary LOSVD tp(ii, e, ui; v p ). In full lines, 
the LOSVD is plotted for u = 90° and e = (black), e = 0.6 
(dark grey), e = 0.8 (grey) and e = 0.9 (light grey). The other 
parameters are set by the choice K = 1. In dashed lines, the 
LOSVD is plotted for ui = 45° and the same eccentricities. 



employed in this context. Hence, we give analytical expres- 
sions for the velocity moments of the binary stars and for 
the LOSVD itself (approximating all binaries to consist of 
stars on circular orbits). 



2 DEFINITIONS 

We first give a number of definitions of important quanti- 
ties and clarify some notations. We focus our attention on a 
single binary system. One star, the "primary", has a mass 
M and the other, the "secondary" has a mass m. All orbital 
parameters pertain to the primary's orbit : 

• v p : the velocity of the primary with respect to the 
binary's center of mass, projected onto the observer's line of 
sight 

• fi p : the gravitational constant of force for the primary's 
orbit around the center of mass. If /j, — G(m + M) is the 
force constant of the orbit relative to the secondary, then 



Hp = 



( 171 V 
\m + Mj 



H = G 



(m + M) 2 



(1) 



• a p : semi-major axis of the primary's orbit with respect 
to the binary's center of mass. If a is the semi-major axis of 
the primary's orbit with respect to the secondary, then 



\m + MJ 



(2) 



.m + M, 

• T, e, i, uj, £1, : respectively the period, eccentricity, in- 
clination, argument of periastron and longitude of ascending 
node of the primary's orbit. The latter does not affect the 
primary's position projected onto the line of sight nor its 
line-of-sight velocity and hence does not enter the calcula- 
tions. 

• 4> '■ the true anomaly along the primary's orbit. 



3 THE SPECIFIC BINARY LOSVD 

The specific binary LOSVD gives the probability of finding 
the primary star with mass M with a line-of-sight velocity in 
the interval [v p — Av p /2, v p + Av p /2], given that it revolves 
around a secondary star with mass m on an elliptical orbit 
with orbital parameters a, e, uj and i. Intuitively, this prob- 
ability is proportional to the fraction of its period during 
which the primary has a velocity in this interval. Mathe- 
matically, this translates into the following expression for 
the specific binary LOSVD : 



(3) 



1 


dvp 


T 


dvp 




dt 



with Vp the line-of-sight velocity of the primary at phase 
angle (f> on its orbit 



Vp = K(ecosuj + cos((/> + uj)). 
Here, 

2na v sin i 



k — 



I /j,p sin i 



Ty/T- 



Hence, if follows that 



dt 



2ttk 



(l_ e 2)3/2 T 

and consequently 

(l_ e 2)3/2 



sin(<^> + uj) \ (1 + e cos ( 



ip(K,e,(j;v p ) = 



2-kk |sin(0 + uj)\ (1 + e cos c 



(4) 



(•») 



(6) 



(7) 



The specific binary LOSVD is a function solely of v p . To 
eliminate the phase angle <j> from this expression, equation 
(^) must be inverted : 



cos(0 + uj) — — — e cos tj. 



(8) 
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Depending on the sign of sin(0 + tj), there are two solutions, 
<p + and if>- : 

1 (l-e 2 ) 3/2 
ip±(K,e,uj;v p ) = , x 



(k) 



2ttk VI- V 2 



^l + eVcoso;±eVl- V 2 



(9) 



with V = v p /k — ecosu). This is a consequence of the fact 
that a star on an elliptical orbit will obtain the projected 
velocity v p twice during each revolution. The LOSVD is then 
just the sum of ip+ and ip- : 

1 (l_ e 2)3/2 

ifi(K,e,Lj;v p ) = , x 



■KK vi - ya 
1 + e 2 sin 2 lu + 2eV cos u> + e 2 V 2 cos(2u;) 
(l - e 2 sin 2 u + 2eV cosw + e 2 V 2 ) 2 



(10) 



In the case of a circular orbit - i.e. e = — this reduces to 
the simple expression 

1 1 



(P(k;v p ) = 



(ii) 



the LOSVD of a harmonic oscillator. The specific binary 
LOSVD is presented for various values of u and e in Figure 
^| The LOSVD is high around the extreme velocities where 
dv p /dt = (i.e. around v p = «(±1 + ecosw)) and around 
the apocenter velocity v p — (e — 1)kcosoj, a direct corollary 
of Kepler's second law. 



4 THE VELOCITY MOMENTS OF THE 
SPECIFIC BINARY LOSVD 

The n th velocity moment of the specific binary LOSVD is 
denned as 



= I Vp ip(Vp) dVp = — J Vp dt 



2tt 



(i-O 



2n3/2 



(cos(0 + LU) + e COS to) 
(1 + e cos <f>) 2 



-#.(12) 



All powers in the integrand can be expanded, yielding 



\n — k , , 
. 6 COS LJ ) X 



^(;)(cos^(sin.)-Mr^( m 2 ) e 



(cos 0) (sin 



(13) 



For the integral to be non-zero, one condition is that k — l = 
2r with r an integer. The other condition is that k + m — 2r 
must be even. Hence, k an m must both be either even 
or odd. This leads to the following expression for the n th 
velocity moment : 



[fc/2] 

E 



fc - 2r 



(coscj) (Sinai) x 



m>0 



-2 



m+n — fc 



r(£±f±i - r)T(r + 1) 



(14) 



Here, [fc/2] stands for the largest integer that is less than 
fc/2 and ^m>o indicates that the summation runs only over 
those m that have the same parity as fc. F is Eu ler's Gamma 
function (see e.g. Gradshteyn & Ryzhik (1965), also for the 
other special functions appearing in this paper). Interchang- 
ing the summations, this expression can be rewritten in a 
more elegant form 



fin 



(k cosw) r 



(1-e 2 ) 



k=0 ' m>0 



(k) 



p / m + fc + 1 N 
m+n— fc \ 2 > 



r(^±2t + i) 

/ fcl-fel-A-m 2 
2Fl l"2'^- ; 2^ ; " tan W 

For e = 0, only the terms with m + n — fc 



(15) 



survive. 

This condition can only be satisfied for even moments and 
if m — and fc = n. Thus, if we substitute 2n for n, and 
make use of the fact that 



2n 

COS U) 2 



F\ ( — n, n; n; — tan 2 o;^ = 1, 

V 2 '2 ' / 

then (hq) reduces to 



1 / ftp . .2 
fj,2n = 3 I — sin % 

1p 



with B Euler's beta function 



(16) 



(17) 



5 THE BINARY LOSVD FOR CIRCULAR 
ORBITS 

The binary LOSVD <p(M ; v p ) gives the probability of finding 
the primary star with a line-of-sight velocity in the interval 
[v p — Avp/2, v p + Avp/2]. To obtain this quantity, the spe- 



cific binary LOSVD fllO] ) must be averaged over all orbital 
parameters. As the calculations for the general case are quite 
cumbersome, we only calculated the binary LOSVD in the 
approximation that the members of all binaries revolve on 
circular orbits. This is clearly not a very realistic assumption 
but it is instructive to have at least this first approximation 
of the LOSVD. To obtain the binary LOSVD, the specific bi- 
nary LOSVD (0) 

must be averaged over all possible values 
of the inclination i, the orbital radius a p and the secondary 
mass m. To do this, we have to assume distributions for each 
of these parameters. 

The orbital planes of all binaries can be expected to be 
distributed randomly. Each direction of the normal of the 
orbital plane is then equally plausible. Using the notations 
of Figure [j], an elementary solid angle around the orbital 
plane's normal is given by dS = sinididQ. Since the direc- 
tion S is uniformly distributed 

1 



47T 



• sin i di dQ. 



<p(S) dS=—dS 
one obtains 

,2* 1 

tfi(i) = / <p(S) dfl — - sini with i £ [0, ty] 
Jo 2 



(18) 



(19) 
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Figure 3. The distribution of the inclination i. The histogram 
shows the distribution extracted from Soederhjelm (1999) with 
the error-bars measuring the Poisson noise. Over-plotted (dotted 
line) is the expected distribution (|19[) . 



3 

87m 



(22) 



with n the average number density of the stars inside the 
radius r. For reasonable values of n, r max is very large. For 
instance, n = 1 pc -3 leads to r max w 10 5 AU. Binaries 
can also be disrupted by an encounter with a third star. For 
typical number densities, stars approach each other to about 
1500 AU once in a Hubble time as they orbit in a galaxy's 
gravitational potential well. Of course, not every such an 
encounter necessarily breaks up a binary pair. Hence, r max 
is at least a couple of thousand AU and can be as large as 
100,000 AU. Our results turn about to be rather insensitive 
to the value of r max . 

If the members of a binary are very close together, one 
of them can be in the other's Roche lobe and be disrupted 
by tidal and centrifugal forces. Back-of-the-envelope calcu- 
lations based on typical masses and radii of K-giants (with 
a mass around 1 Mq and a radius of approximately 10 R©)- 
which one is most likely to observe since they are so luminous 
- give r m i n m 15 R© as a fair mean value. 

Due to 



as the distribution of the inclination. In Soederhjelm (1999) 



<p(a P ) da p 



(23) 



the orbital parameters and masses of about 200 nearby vi- 
sual binaries are given. The data are obtained from Hippar- 
cos and ground-based telescopes. Figure ^| shows the distri- 
bution of the inclination, extracted from this catalog, and 
clearly demonstrates the validity of (|l9|). 

The distribution of the semi-major axis cannot read- 
ily be extracted from a catalog. Catalogs of visual binaries 
usually only contain the semi-major axis expressed in arcsec- 
onds since their distances are not known. Catalogs of spec- 
troscopic binaries do give the semi-major axis expressed in 
e.g. AU but they tend to be strongly biased towards very 
tight binaries since these will have the highest orbital ve- 
locities and hence the largest and the easiest measurable 
Doppler-shifts. However, it is natural to assume that more 
tightly bound binaries are less easily disrupted by external 
influences. Hence, one would expect to observe less binaries 
with lower binding energies. We therefore adopt 



<p(a) da 



E~ 



' da 



(20) 



as the distribution of a (the orbital radius with respect to 
the secondary). Here, E is the binary's binding energy and 7 
is a real number. The closest gravitational few-body system 
at hand, our planetary system, obeys (|2(]) quite well with 
7 «1. 

The upper and lower bounds for a, r m i n and r max re- 
spectively, can be estimated as follows. For a binary to be 
stable, the mutual gravitational attraction between its mem- 
bers should overcome the tidal forces exerted by the host 
galaxy. In other words, if we denote the average mass of a 
star by m : 



n ■ 



2M(r) 



(21) 



with r the distance of the binary to the galaxy's center and 
M(r) the total mass inside a sphere with radius r. This leads 
to the following expression for the maximum orbital radius : 



with a m = m/(m + M)r m i n and analogously for om. 

For the distribution of the secondary mass, we take the 
Salpeter initial mass function 

1 - x 



tp(m) dm - 



dm 



(24) 



with x = 2.35. The lower bound for the mass distribution 
is taken to be mi = 0.08 Mq, the mass of the most light- 
weight stars capable of nuclear fusion. The upper bound m2 
is treated as a free parameter of the binary orbital distri- 
bution. Realistic values for mi are around 1.25 M©. This 
corresponds to stars with a total life-time of 5-6 Gyr, com- 
parable to the ages of the stellar populati ons in Local Group 
dwarf galaxies (e.g. Smecker-Hane et al. (1994)). These dis- 
tributions can be c ompared to those found by Duquen- 
noy & Mayor (1991) for a sample of 164 solar-type stars 
in the solar neighborhood. For their sample, they find the 
following distribution of the binary periods : 



<^(log T) dlog T oc exp 



log T - log T 

flog T 



dlogT(25) 



with logT = 4.8 and 0"i og T = 2.3 (T is expressed in days). 
Periods were found to lie in the range — 1 < log T < 10. Us- 
ing @, @ and r max w 100, 000 AU (applicable if binaries 
are disrupted by tidal forces exerted by the host galaxy) , we 
find an approximately uniform distribution of logT in the 
range < log T < 10. Using r max « 2000 AU (applicable if 
binaries are disrupted by a close encounter with a third star), 
yields an approximately uniform distribution of log T in the 
range < logT < 7.5. As always, the truth lies somewhere 
in between these two extremes. The Salpeter law adopted 
here provides a reasonable fit to the mass-ratio distribution 
although it over-estimates the number of low-mass compan- 
ions with q < 0.4 found by Duquennoy & Mayor. 

We start with the integration of the specific binary 
LOSVD (|ll|) over the inclination, making use of the dis- 
tribution (191). The integration region must be limited to 
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arcsin to < i < 7T — arcsin to 



(26) 



with to 



-£-£. to assure that the argument of the square 

root in the denominator of is positive. This, of course, 
is because, for a given a p and \i p , some binaries with a large 
inclination cannot reach the velocity v p and must be ex- 
cluded from the integration. The binary LOSVD after this 
integration becomes 



<p{M-v p ) = 



2lYV v 



7r — arcsin £q 



1 ■ 2 • 

-5 sin i 



■. di 



(27) 



We employ the distribution function ( |23[ ) to average for- 
mula ( p7[ ) over the semi-major axis. The integration region 
is now limited to 



a m < cip < mm am 



/')--■ 



(28) 



This ensures the positivity of the binary LOSVD by exclud- 
ing all binaries that, for a given fi p , cannot reach the velocity 
tip. The binary LOSVD at this stage of the integrations takes 
the form 



<fi(M;v p ) = 



1 



'7 



27 r 



1-7 

max ' 



mm am 



Hp 



(m + A/) 2 ' x 
(29) 



The condition am < ^ is satisfied if 



m > mo = 



GM 



+ 



''max 

GM 



4 + 



GM 



(30) 



In that case, there is no danger of the integrand becoming 
negative. If m < mo however, the integration region must 
be limited to 



M 

m > 7713 = — 



GM 



+ 



GM 



GM 



(31) 



if 7 < 3/2. Clearly, 7713 << mo if v p > 0. The integration 
over the secondary mass involves three basic integrals. The 
first one is 



h(xi,x 2 ) = 



•7 



3/2-7 



3-2 7 i£2-riT VGml 



x 1 M 



1 



•7 



3/2-7 
' max 



L (m + M) 1/2 dm 

1 1 - IE 



3-2 7r 



1-7 



rillVGml- 



M 



1/2- 



' /•' ! —x — 1, ^;xi,X2 



(32) 



with 

F(a, 6;a;i,a;2) 



ir a (l + x) 6 dx, IEl < X2 

(a, b; x\, X2) if a;2 < 1 
.F(a,&;a:i, 1) + 



£2 

if 2:2 > 1 but x\ < 1 

= ^(-a-6-2,6;— , — 
£2 £1 , 

if X2 > 1 and iei > 1 (33) 



■(- 



and 



F (a,b;x 1 ,x 2 ) 



a+1 



>Fi(-b,l + a;2 + a;-x 2 ) 



—^TJ 2^i(-fc, 1 + a; 2 + a; -si). (34) 

Here, 2-F1 is the hypergeometric function. The second inte- 
gral is completely analogous : 



h(xi,x 2 ) 



1 — 7 r 



3/2- 7 



1 



1 -x 



max 



m — m. 



m-"- 1 (m + M) 1/2 dm 



, 3/2-7 , , 

1— <y 1 — 7 /T^y 1 — x Im 1 — x 



3-2 7r 



■ n 7 VG m\- x - m} 



M 1/2 -"J , (-i-l,i;xi,aai 



(35) 



The third basic integral is 



h(v p , xi,x 2 ) = 
J_(G\ 

Vg W p 



■7 



3 - 2t r 1-7 - 

/ ' max ' m j 



1 — 7 ~vi 1— 



3/2-7 /.zoM 



-27 — x+2 



1-7 



(m + My 1 dm 
1-x 



1— -y 1 — 7 ™1— a; 
" - r min m 2 



3 - 2 7 r n 

M x/a— / GM \ 3/2-7 

!f££ F(-2 7 -x + 2,7-l;x 1 ,x 2 ). (36) 

VOr \ V P J 

In each of these integrals, the factor 

1-7 1 M 1 / 2 -* I-3; 



Q 



3-27 r ^ a 7 - r 1 ; 7 \ < r - , 



s/G m\- x -m 1 - :l 



(37) 



can be taken out. We are then left with the following three 
building blocks : 



1i (xi,x 2 ) 
l 2 (xi,x 2 ) 

Ja(v p ,xi,X2) 



.3/2— y 



F l-x 



l, -^\xi,x 2 



3/2-7 „ / , 1 

- r Jnn 1 F\-X-l,-; Xl ,X 2 

F (-27 -a; + 2, 7-I; 3:1,3:2) 



(38) 



The binary LOSVD, averaged over all orbital parameters, 
finally takes the following form : 



<p(M;Vp) 



Q 



X\ ( max 



J mi m l 

l¥' Uj 



X2 I max ■ 



m 2 
U 

+ 



+ 



f m 3 mil m 2 
I M ' ¥/ ' M , 
I ( m-z mi "I fm2 mol\ 

(^,max|— ,— j,mm(— ,— jj 



(39) 



As can be expected, the LOSVD goes to zero for m$ > m2 
or equivalently for 
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Figure 4. The logarithm of the binary LOSVD <fi(M; v p ) for 
circular orbits. The primary mass is M = 0.8 Mq, the secondary 
mass is averaged over the interval 0.08 — 1.5 Mq, r m i n = 15 Rq, 
r max = 101559.1 AU as 0.5 pc (n = 1 pc" 3 ), 7 = 1 and x = 
2.35. The inset shows the flat top of the LOSVD for very small 
velocities. For the parameter values adopted here, the LOSVD is 
constant for velocities smaller than »i fs 7 m/s. 



Vp > ^3 = 



Gm\ 



(40) 



(m 2 + M)r min ' 

the highest line-of-sight velocity any binary can reach. As 
long as mo < mi, or equivalently for 



v„ < Vi 



(mi + M)r r 



(41) 



the LOSVD is velocity independent. This is the highest ve- 
locity that the slowest binaries, i.e. the ones with the lowest 
secondary mass and the largest orbital radius, can reach. 
Usually, r max is very large and consequently, vi will be a 
very small velocity, of the order of a few m/s (see Figure^). 

To be exact, the binary LOSVD ( p9[ ) still needs to be 
multiplied with the probability that the primary star is actu- 
ally observed and then be averaged over the primary mass. It 
is this final function that quantifies the extra broadening of 
the absorption lines due to the binary population. However, 
the probability of observing a star of a particular spectral 
type depends both on its intrinsic brightness and on the 
number of such stars that are around. As the brightest stars 
are also the most scarce ones while the faintest stars are ex- 
tremely abundant, this probability will be sharply peaked. 
For old to intermediate stellar populations, this peak corre- 
sponds to early K giants on the tip of the giant branch. We 
will therefore not do this last integration and instead use 
an average value for the primary mass. Most authors use 
M = 0.8 Mq and this is the value we will adopt here. 



6 THE MOMENTS OF THE BINARY LOSVD 
FOR NON-CIRCULAR ORBITS 

The binary LOSVD for elliptic orbits is not analytically 
tractable, but its velocity moments, which we denote by Jin, 



0.2 0.4 0.6 0.B 1 

eccentricity e 

Figure 5. The distribution of the eccentricity e. The histogram 
shows the distribution extracted from Soederhjelm (1999) with 
the errorbars measuring the Poisson noise. Over-plotted in dotted 
lines is the distribution (^j). The uniform distribution is plotted 
in dashed lines. 



are. We set out with the expression ( [lil ) for the moments of 
the specific binary LOSVD and average it over the argument 
of the periastron u), the inclination i, the secondary mass m, 
the semi-major axis a p and the eccentricity e. For i and m, 
one can use the distributions ( |l9|) and (pi|). The distribution 
( p3| ) for the semi-major axis is also still applicable but care 
has to be taken with the integration limits : the bounds for 
a p are set by the conditions 



a(l + e) < Tmax, 

a(l - e) > r min , 



(42) 



which makes the integration limits of a p dependent of e. 

As can be seen in Figure the distribution of the ec- 
centricity of a sample of 200 visual binaries can be approxi- 
mated fairly accurately by 

(1 + a)e a (l - e 2 ) 13 

~ v ' a 4-1 7-1/ n 1+rv 9 \ V / 



2 



3+q . 2 

9 ; '-max 



with a = 0.5, P = 1 and maximum eccentricity e max = L 
The e ccentricity distribution found by Duquennoy & Mayor 
(1991) shows that tight binaries have more circular orbits 
than wide ones but, overall, their eccentricity distribution is 
similar to the one derived by Soederhjelm. In the remainder, 
we consider the parameters a, (3 and e ma x to be fixed by the 
observations. This distribution reaches a maximum for 



('to 



2/3- 



0.45, 



(44) 



for the parameter values adopted here. The uniform distri- 
bution 95(e) = 1/emaxj, which is the distribution advocated 



by Hargreaves et al. ( 1996 ) and Mateo et al. (1992), is ob- 
tained by setting a = P = 0. 

The argument of the periastron u) is obviously dis- 
tributed randomly between uj — and uj — 2n. Therefore, 



<p(iu) duj = — — dui. 
2n 



(45) 



Averaging the n* 1 velocity moment of the specific bi- 
nary LOSVD (|l5|) over us involves the integral 
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l / 

— / (cosoj)™ - r (sinu;) r du> (46) 
2n Jo 

which can only be non-zero if n is even. We will therefore 
consider only the even moments and make the transforma- 
tion n — > 2n in which case 



2tt 



(COS LU) (SinU) did — 

1 - r + l/2)r(r + 1/2) 

7T 



(47) 



T(n + 1) 

The integrations over i and m are pretty straightforward 
and one obtains the following expression for the moments 
fl2n ■ 

1 - x 



H2n = 



/ GM \ n 1_ 



V a / (2. + l)0F (ff ) 1 -_(^) 1 



F ^2n — cc, — n, 



2n / s (fe) 
JTll "i2 \ \ I 2n\ \ ■* I —2 



Ei E 



MM) \ k I * — ' \ m 

fc=0 x ' m>0 



m + 2n— fc 



(i-O 



2n3/2-k 1 V 2 , 



}F 2 



' k + m 
~ 2 

fc 1-k 1 1-fc-m 1 



r /«±m +1 ] 



2' 2 '2' 



n: 1 



(48) 



The next step is the integration over the semi-major 
axis a, making use of ( pfj| ) : 



1-7 



,1-7 _ „l-7 



1 + e 1 — -7 

a da — 



' ' max ' m i n 



1 — 7 — n 



1 - e 



r," 

1 — 7 — n 



(49) 



Finally, we are left with the integration over all possible 
values of the eccentricity. Making use of the lemma 



Q(a,b,c-x) 



y a (l+y) b {l - yfdy 



a + 1 



(-b) m (a+l)r . 

(a + 2) m m! 



2-Fi(— c, a + m + 1; a + m + 2; x) 

one finds that : 

(GM) n (l-j)(l-x) 1 



(50) 



A*2n 



R (l-7-n)(2n + l) 



1 — 7 1 — 7 



1 + a 



\~W) ~\MJ m " 2 H P ' 2 ' 2 ' e maxj 

m 
"'Ik 



, mi m 2 \ 

F|2n-z,-n;-,-]x 



fc = x ' m>0 



r(*4p + i) 



(4 



, fc 1 — fell — fc — ml 
3-F2 I --, — — , -; , - - n; 1 ] x 



2 ' 2' 2 '2 



^1 — 7 — n 



C/(m + 2n — fc + a, - +7 + /?, - — n + ,3;e n 



1 — r-n. 



r min ( m + 2n - fc + a, - - n + /?, 



- + 7 + /3;e n 



(51) 



For e ma x = 0, only the terms with m + 2n — fc = survive. 
This is only possible if m = and k — 2n and ( |5l| ) reduces 
to 



^2n = 



(GM) r 



■ 7 



„l-7-n 



^1 — 7 — n 



2n 



1 1 — 7 — n 

, x ~ x , f( 



r 1 — ' _ r 1 ^ 7 
' max ' min 



2n-x,-n;^,^), (52) 



the expression for the moments of the LOSVD for circular 
orbits. 

In the next paragraph, we discuss how the velocity dis- 
persion of the binary LOSVD for non-circular orbits depends 
on the distributions of the different orbital parameters. An- 
other important point that will be addressed is the way the 
shape of the observed LOSVD depends on the binary frac- 
tion. 



7 DISCUSSION 

7.1 The binary velocity dispersion as a function 
of the orbital parameters 

As is obvious from Figures |^ to [To| the binary velocity dis- 
persion, i.e. the second order velocity moment fi2, depends 
most sensitively on the inner cutoff radius r m i n and the ex- 
ponent 7 in the distribution of the semi-major axis a v and 
to a lesser degree on the primary mass M and the exponent 
x in the Salpeter IMF adopted as the distribution of the sec- 
ondary mass m.2 . As the stars in a binary system are allowed 
to approach each other more closely - in other words, as r m i n 
is lowered - their orbital velocities rise rapidly, boosting the 
velocity dispersion of the binary population as a whole to 
values as high as 10 km/s for r m i n = 1 Rq, as can be seen in 
Figure ^[ Since the bright giant stars that one is most likely 
to observe have radii of the order of 10 Rq, the velocity dis- 
persion is limited to values below 4 km/s. The value of the 
outer cut-off radius, r max , does not strongly affect the value 
of the binary dispersion. Higher r max values mean that more 
stars orbit on wide and consequently slow orbits, causing the 
dispersion to be a bit lower. A higher value of 7 causes stars 
to move on more tightly bound orbits, thus boosting their 
velocities. This is clear from Figure Nfl. If x, the exponent 
of the Salpeter IMF, is made larger, the secondary masses 
will tend to be lower which also causes the primary's orbital 
velocity and hence the velocity dispersion of the binary pop- 
ulation as a whole to be lower (see Figure [ll]). If M is raised, 
the center of mass will be closer to the primary, causing its 
orbital velocity to drop slowly, as is observed in Figure ^| 
Larger values for m2 will produce heavier secondary stars. 
A slowly rising velocity dispersion ensues, as can be seen in 
Figure H. The same exercise was done with a uniform dis- 
tribution for the eccentricity. The outcome was essentially 
identical : the binary velocity dispersion differed at most by 
0.1 km/s. Hence, our results do not depend critically on the 
adopted distribution for the eccentricity. 

The values for the binary velocity dispersion obtained 
here are about the same as those obtained by Hargreaves et 
al. although these authors use Gaussian distributions for the 
ellipticity, secondary mass and period, making it difficult to 
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Figure 6. The binary velocity dispersion as a function of the in- 
ner cutoff radius, r m ; n . The values of the parameters are indicated 
in the plot. 
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Figure 7. The binary velocity dispersion as a function of the 
outer cutoff radius, r max . The values of the parameters are indi- 
cated in the plot. 



directly compare the results. Mateo et al. ( 1993] ), employing 
a Monte Carlo simulation with a uniformly distributed ellip- 
ticity and secondary mass and a power law distribution for 
the period, find higher velocity dispersions. However, it is 
mentioned by Hargreaves et al. and by Olszewski et al. that 
the values published by Mateo et al. are overestimated due 
to a coding error. The corrected Mateo et al. values however 
are in good agreement with the results obtained by other 
authors and with those presented here. 

In the following, we will adopt the model with 



7 
x 

^min 

m i 

1712 

M 



1 

2.35 
15 R Q 
0.08 M 
1.25 M 
0.8 M 



(53) 



as our "standard" model. Its velocity dispersion is Ob = 
2.87 km/s and its kurtosis amounts to £4^ = 86.89. 




m 2 =1.25M e 



r min- l5R 



0.5 



M (M ) 



Figure 8. The binary velocity dispersion as a function of the 
primary mass M . 



7.2 The observed LOSVD 

The LOSVD that would actually be observed is the convolu- 
tion of the binary LOSVD (describing the motion of stars in 
binary systems) and the intrinsic LOSVD of the galaxy (de- 
termined by the stars' orbital motions through the galaxy). 
If we denote the intrinsic galaxy LOSVD by fi(v p ) and the 
binary LOSVD by ipb(v p ), then the observed LOSVD tp {v p ) 
can be written as 



<fii(v p — x) ((1 — a) 5{x) + a ifib(x)) dx 



(1 — a) <fii(v p ) + a / <pi(v p — x) tfib(x) dx (54) 



with S(x) the Dirac delta function and a the binary frac- 
tion, i.e. the fraction of visible (solar-mass) giant stars that 
have a companion. As an indication, Duquennoy & Mayor 



(1991) find a « 0.6 for solar- mass stars in the solar neigh- 
borhood with the same range of binary periods as adopted 
here. In Figure [t], the observed LOSVD is plotted for dif- 
ferent binary fractions. The intrinsic LOSVD is a Gaussian 
with dispersion cti = 3km/s. The binary LOSVD is given by 
@ with M = 0.8 M , r min = 15 R , r max = 101559.1 AU 
(n = 1 pc -3 ), 7 = 1 and x = 2.35. The secondary mass 
is averaged over the interval 0.08 - 1.5 M . The LOSVD 
becomes clearly more peaked and develops broader high- 
velocity wings as the binary fraction increases. A compre- 
hensive way of assessing the change in the shape of the ob- 
served LOSVD as a function of the binary fraction is to 
study its moments. If the intrinsic galaxy LOSVD is a sym- 
metric function of the line-of-sight velocity, then one finds 
for the velocity dispersion a Q and the fourth order moment 
^4 j0 of the observed LOSVD : 



2 . 2 
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Figure 9. The binary velocity dispersion as a function of the 
upper mass cutoff m2. 
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Figure 11. The binary velocity dispersion as a function of the 
exponent x in the Salpeter law for the secondary mass distribu- 
tion. 
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Figure 10. The binary velocity dispersion as a function of the 
exponent 7 in the distribution of the major axis a p . 



H4,o = /U4,i + 01(14, b + 60(7; a b 



(55) 



with o"; and (i4,i respectively the second and fourth order 
moments of the intrinsic galaxy LOSVD and u\ and fj,^ 
respectively the second and fourth order moments of the 
binary LOSVD. The fact that we have the fourth order mo- 
ment of the binary LOSVD at our disposal allows us to 
study not only the observed velocity dispersion but also the 
observed LOSVD 's kurtosis 



?4,o = — 

at 



(56) 



as a function of the binary fraction. 

From Figures [D| ^ and ^| it is clear that the ve- 
locity dispersion of only the most light-weight galaxies will 
be appreciably affected by their binary population. There, 
the velocity dispersion of the observed LOSVD is presented 
as a function of the binary fraction for intrinsic Gaussian 
LOSVDs with <Ti = 3 km/s, o\ = 6 km/s and o\ = 9 km/s. 




Figure 12. The observed LOSVD for different binary fractions 
a, indicated in the plot. The intrinsic LOSVD is a Gaussian with 
dispersion <t; = 3 km/s, the binary LOSVD is given by ( |39[ ) with 
M = 0.8 M Q , r min = 15 R Q , r max = 101559.1 AU (n = 1 pc" 3 ), 
7 = 1 and x = 2.35. 



The kinematics of the binary population are those of the 
"standard" model. Only if the intrinsic velocity dispersion 
is of the same order as the binary velocity dispersion, i.e. 
<Ti ~ 3 km/s, is there a noticeable change in the velocity 
dispersion. A binary fraction of 60% raises the velocity dis- 
persion of the cr; = 3 km/s LOSVD with 24%, that of the 
cr; = 6 km/s LOSVD with 7% and that of the a\ = 9 km/s 
LOSVD with 3%. These Figures, however, yield a surpris- 
ing result : whereas the influence of the binary population on 
the observed velocity dispersion is limited to stellar systems 
with intrinsic dispersions that are comparable to the binary 
velocity dispersion, the shape of the observed LOSVD - as 
quantified by the kurtosis - is appreciably affected for sys- 
tems with intrinsic dispersions as high as three times that 
of the binaries. The kurtosis can be related to a quantity 
that is of more direct interest to observers : the coefficient 
of the fourth order term in the expansion of the LOSVD 
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Figure 13. The velocity dispersion and kurtosis of the observed 
LOSVD as a function of the binary fraction a. The intrinsic 
galaxy LOSVDs is a Gaussian with cr; = 3 km/s. The parameter 
values are those of the "standard" model, i.e. with <r b ss 3 km/s. 
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Figure 15. The velocity dispersion and kurtosis of the observed 
LOSVD as a function of the binary fraction a. The intrinsic 
galaxy LOSVDs is a Gaussian with <jj = 9 km/s. The parameter 
values are those of the "standard" model, i.e. with <r b ss 3 km/s. 
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Figure 14. The velocity dispersion and kurtosis of the observed 
LOSVD as a function of the binary fraction a. The intrinsic 
galaxy LOSVDs is a Gaussian with Ci = 6 km/s. The parameter 
values are those of the "standard" model, i.e. with <rt, ss 3 km/s. 



in Gauss-Hermit e pol ynomials, hi (Gerhard ( 1993 ), van der 
Marel & Franx ( |l993| )) : 



hi 



e 



4,o - 



8^6 



(57) 



E.g. a binary fraction a = 0.6 raises the kurtosis of a o\ = 
6 km/s Gaussian to 5.09 (hi ss 0.1) and that of a <7i = 9 km/s 
Gaussian to 3.47 (hi w 0.025). For stellar systems with lower 
intrinsic velocity dispersions, the LOSVDs will be even more 
distinctly non-Gaussian. 



8 CONCLUSIONS 

The maximum additional velocity dispersion due to binary 
stars is estimated at <Tb w 3 km/s, in good agreement with 
other authors. Only stellar systems with intrinsic velocity 
dispersions comparable to this value will have observed ve- 
locity dispersions that are noticeably affected by the pres- 
ence of binaries, i.e. dwarf Sphero idals such as those found 
in the Local Group (Dekel & Silk (1986)), globu lar clu sters, 
low-surface-brightness disk galaxies (Bottema (1993)) and 
the central regions of nucleated dwarf ellipticals. Not only 
the velocity dispersion but also the shape of the galaxy's 
LOSVD is altered. The presence of binaries has a measur- 
able effect on the kurtosis of the observed LOSVD for stellar 
systems with a velocity distribution as high as 10 km/s. For 
stellar systems with lower velocity dispersions, the LOSVDs 
will be even more distinctly non-Gaussian, being strongly 
peaked with broad wings. This feature can mimic radial 
anisotropy which combined with the enhanced velocity dis- 
persion could lead to an over-estimation of the dark matter 
content of dSphs and globulars. 

If the kinematics of a galaxy are derived from the radial 
velocities of discrete stars - as is the case for the dwarf 
spheroidals in the Local Group - repeated observations can 
eliminate the binaries from the star sample. Methods that 
rely on integrated light spectra - e.g. the crowded regions of 
globular clusters or dwarf galaxies outside the Local Group 
- will be plagued by the above mentioned effects. 
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